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Abstract. We investigate experimentally the distribution of the free volume inside 
a bidimensional granular packing associated either to single grains or to clusters of 
grains. This is done for two different kinds of grains and two levels of compaction. 
The logarithm of the free volume distribution scales in a nonextensive way with 
the cluster size for cluster sizes up to a few hundred grains. Having factorized this 
size dependence, we characterize the distributions by two intensive parameters. We 
discuss the interpretation of these parameters and their possible relation to Edwards' 
compactivity. 
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1. Introduction 

Granular media consist in a large number of grains and one would like to have 
general predictions for the macroscopic laws ruling such media without having to deal 
with individual components. Recent experiments suggest that such a thermodynamic 
description might be feasible |2], and there is a strong motivation for providing a 
statistical ground to this hypothetic description [3J EJ 13 El HI IE] • Yet, a theoretical 
understanding is still lacking regarding the way microscopic configurations are explored 
dynamically. 

It has been argued by Edwards and collaborators |HJ IHj that the dynamics is 
controlled by the mechanically stable -the so-called blocked- configurations and that 
all such configurations of a given volume are statistically equivalent provided that the 
driving involves extensive manipulations, such as shaking, shearing or pouring. This 
immediately leads to the definition of a configurational entropy S(V) which counts the 
number of such configurations compatible with a given volume. The associated state 
variable, the "compactivity" X, then naturally follows from the relation X' 1 = dS/dV. 
So far, this equiprobability assumption has been tested through numerical simulations 
in different physical contexts: in schematic finite-dimensional models with kinetic 
constraints @J, in spin glass models with non-thermal driving between the blocked 
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states ^U] and finally in a few more realistic models of particle deposition jTTj or 
molecular dynamics simulations of granular media driven by an external shear [7j . Note 
however that Edwards' approach has also been shown to fail in some cases |12[ llrij. 

^From an experimental point of view, a stationary dynamics has been observed 
for vertically vibrated glass beads at rather high volume fraction, and has been shown 
to be independent of the history of the system pQ. In addition, a relation between 
the fluctuations of volume fraction and the volume fraction itself has recently been 
reported =2J in a similar experimental system. Assuming that usual thermodynamic 
relations are valid, this led to a determination of the compactivity X of the packing. 
An alternative determination of X has also been proposed from the exponential decay 
of single grain volume distribution ^] . 

In spite of these globally encouraging results, the proposal made by Edwards and 
coworkers still needs to be clarified in several respects. First, the choice of the physical 
variables to be used in practical situations essentially remains an open issue JHI, and 
the computation of the corresponding density of state turns out to be out of reach of 
present theories [T7j. Second, this statistical construction may be strongly affected by 
the presence of correlations that extend over length scales comparable to the system 
size. Indeed the occurrence of such long-range correlations has been reported in various 
contexts, for instance in jammed granular systems [T%1 in porous media or in 
cellular structures j2I] • In addition, recent theoretical advances suggest that amorphous 
systems are close to a critical point associated with a jamming transition [22112312115 thus 
also supporting the existence of correlations on length scales small than a correlation 
length that may be rather large. These various issues clearly call for experimental checks 
of the foundation of Edwards' approach. 

In this paper, we report on an experimental investigation of the distributions 
of volume, as well as free volume, of Voronoi cells associated with each grains in a 
bidimensional granular packing. The distribution P(v{^) of the free volume per grain 
v N in clusters of neighboring grains is also investigated in a systematic way as a function 
of the size N of the clusters. Two different types of grains, as well as two different levels 
of compaction, are successively used. Our central result is that lnP(f^) scales in a 
nonextensive way as N a g(v^ [ ,Xe,Ve) in the accessible range of N, with a ~ 0.8. The 
parameter Xe and r) e are intensive, and Xe niay be thought of as a generalization of the 
compactivity introduced by Edwards. 

The paper is organized as follows. In Sect. 2, we briefly recast Edwards' formalism 
into a more general form which allows possible departures from the equiprobability of 
states to be taken into account. The experimental set-up is presented in Sect. 3, and 
the experimental results on the distributions of volume and free volume are reported 
in Sect. 4. Finally, Sect. 5 is devoted to a discussion of the physical interpretation and 
consequences of our results. 
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2. General theoretical formulation 

In order to propose an experimental procedure for the determination of statistical 
properties of granular packings, we first discuss in more details the theoretical 
background. Although the definition of compactivity apparently relies on the 
assumption that the states compatible with a given volume are equiprobable, a 
compactivity can actually be defined in a more general context. Let us consider a 
granular packing, and denote its microscopic configurations by C. In the spirit of 
Edwards, we assume that the volume is a globally conserved quantity. On general 
grounds, the steady-state probability P{C) may be written as, taking into account the 
conservation of volume 

P(C) = ^1 ) 5(V(C)-V) (1) 

which defines the probability weight f(C) and the normalization factor Z(V). Let 
us assume that, on large scale, different regions of the system are only correlated by 
the global volume conservation. By this we mean that, considering a partition of the 
system into two subsystems S± } 2, Z(V) factorizes into Z(V) ~ Zi(Vx) Z 2 (V — Vi) -the 
approximate equality results from possible interface contributions, that are expected 
to vanish in the thermodynamic limit. This is in particular the case if the probability 
weights f{C) are assumed to factorize as f{C) = /i(Ci)/2(C2), where C = (Ci,^). 
Under such assumptions, a compactivity can be defined as a generalized "intensive 
thermodynamic parameter" for out-of-equilibrium systems |25j : 
1 d\nZ 

X ~ dV ' { ) 

Note that Edwards' definition is recovered when f(C) is uniform over blocked 
configuration, and zero for other configurations. This reformulation shows that the 
equiprobability assumption is actually unnecessary to define a compactivity. 

A relevant experimental situation consists in studying a small part of a large 
packing. In this case, the effect of the rest of the system, acting as a reservoir of volume, 
may be accounted for by its compactivity. Integrating out the degrees of freedom of the 
reservoir in distribution ((TJ) yields the "canonical" distribution for the subsystem 

p ^) = Y§cf v{c)lx (3) 

Interestingly, due to the presence of the exponential factor in Eq. Q, usual 
thermodynamic relations like 

{5V2) = (4) 
are still valid independently of the specific form of the probability weights f(C), provided 
that the above factorization assumption is satisfied. Such thermodynamic relations offer 
the possibility to measure the compactivity of a given packing, as done in |2 . Yet, more 
information may be gained by measuring the probability of observing a volume V in a 
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subsystem of N grains. This distribution <&n{V) is obtained by summing the canonical 
distribution Q over all configurations with a given volume V: 

*n(V) = J dCP c (C)S(v(C)-v) = Y^) e ~ V ' X = zJx) eSNiV) ~ V/X{5) 

Apart from the exponential factor, most of the information about the system lies in the 
function Zn(V), the logarithm of which is a generalized entropy Sn{V) = lnZjv(V). 
Computing Sn(V), which depends in a complicated way on the (usually unknown) 
probability weights f(C) is a formidable task, out of reach of present theories. Hence, it 
is of primary importance to gather quantitative informations on Sn(V) and its scaling 
in N from the experimental side. 

To this aim, one would first need to choose a set of physical variables that are both 
experimentally accessible, and relevant to describe a microscopic configuration. These 
variables may actually not fully define the configuration, as different configurations 
might lead to the same values of the variables. The most simple choice is to consider 
the volumes Vi of cells associated to each grain through some tiling of space , leaving 
aside from the description the forces between grains. Alternatively, the free volumes v{, 
defined by subtracting to u$ the volume of the grain, may also be a relevant description, 
more convenient for polydisperse systems. In this paper, we consider both descriptions, 
with still an emphasis on the free volume description. In addition, we focus on systems, 
like our experimental set-up presented below, that are constrained to evolve among 
blocked states. Such systems may then be described by a density of state p({v{}), which 
corresponds to the number of blocked states with given free volumes v{ . Introducing 
the total free volume of the packing V? = V — V g , with V g the total volume occupied 
by the grains, the function Z N {V^) reads 

Z N (V f ) = I ]\dv!p({v{})f{{v{})6 ^ 

Let us comment briefly on the physical interpretation of the above equation. First, 
we note that p({v{}) is the static density of states, which corresponds to a simple 
"counting" of blocked states compatible with the constraints. p({v{}) encodes a 
compatibility condition regarding the variables that are not taken into account in 
the chosen microscopic description - for example the forces - and hence measures the 
degeneracy of the corresponding configuration. Second, the dynamical population of 
the states is accounted for by the dynamical probability weight f({v{}). This factor 
precisely encodes the possible departure from equiprobability. Thus these two terms 
clearly have a different interpretation. Yet, only their product appears in equation (JHJ), 
so that it is not possible to determine f({v{}) without an independent knowledge 
of p{{v{}). As such a determination of p{{v{}) is in general out of reach in a 
realistic experimental system, it follows that the validity of Edwards' assumption on 
the equiprobability of states cannot be tested experimentally. Only in some simplified 
theoretical models p({v{}) may be determined, and the validity of Edwards' assumption 
may be tested explicitly [HE1II3!- Yet, as explained above, Edwards' assumption is not 
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Figure 1. (a) Experimental set-up. (b) Sketch of the modified Voronoi tessellation 
for bidisperse systems. 

necessary to implement the above statistical formulation and to motivate experimental 
determination of the (generalized) entropy and of the associated compactivity. 

3. Experimental set-up 

Let us now describe our experimental set-up, which is presented schematically on 
figure l^a). It is made of a bidimensional rectangular glass container of width 50 cm and 
height 60 cm, with an inner gap of thickness 3.5 mm. This container is filled with a single 
layer of 5000 cylindrical brass spacers of two different sizes in equal number, resulting in 
a homogeneous and disordered bidimensional packing. The spacers (that we shall also 
call "grains" ) occupy approximately half of the total volume of the container, as shown 
on figure [Ha). The two kinds of spacers have a thickness of 3 mm, and differ by their 
diameter, which is respectively d s — A mm for the small spacers, and di — 5 mm for the 
large ones. For commodity, we choose d s as the unit length, so that the respective areas 
of small and large disks are v° s = it/A ~ 0.78 and v® = it/A (di/d s ) 2 ~ 1.23. 

In addition to the two different kinds of grains which differ by their size, we also 
use two different families of grains (each with small and large grains) with very different 
roughness properties. The first family consists in metal grains with a smooth surface, 
the static friction coefficient of which is approximately fi ~ 0.35, whereas the second one 
is composed of serrated rollers with a crenel shape on their peripheral surface, resulting 
in a static friction coefficient /i ~ 0.5. Moreover different densities of the packing are 
obtained by vibrating it with a hammer-like device installed below the container. 

The aim of the experiment is to generate a large number of different random 
configurations of the packing, and to analyze these configurations through an image 
processing. The protocole is chosen in order to achieve a maximal randomness in the 
exploration of all the accessible configurations even though it is difficult to have a proper 
characterization of such a randomness. However we could verify that the individual 
distribution of voids and their repartition in space is very different from one image to 
another. An experimental run consists in 5000 cycles defined by the following procedure. 
Starting from a vertical position, the cell is slowly rotated of 360°, at a constant angular 
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velocity of 6°/s. During such a cycle, the grains fall from one side to the other and 
then back to the initial side. We let the packing relax for a few seconds so that it reach 
a mechanically stable state and we take a picture of a suitably defined window within 
the bulk of the packing, which contains approximately 1000 grains. Once this is done, 
the whole system is vibrated for 10 s at a frequency of 10 Hz and with an amplitude 
of 5 mm, to allow for rearrangements and compaction. A second picture is then taken. 
Note that vibrating the system does not aim at studying the compaction dynamics itself 
but simply at generating a second packing state with a different density. 

The image processing consists in the following steps. The first step is to locate in 
each picture the centers of the grains. Then a Laguerre tessellation is used to associate 
to each grain a specific area. This tessellation is a modified version of the Voronoi one 
with respect to standard definitions, in order to take into account the bidispersity of 
the grains |2Z|. An example of such a tessellation is shown on figure H^b). Once all cells 
are defined, the areas of the cells are computed and recorded, keeping also track of the 
position and type (small or large) of the associated grains. 

Using the above procedure, we obtained three different sets of data, corresponding 
respectively to smooth grains without vibration, smooth grains with vibration, and 
rough grains without vibration. 

4. Experimental results 

4-1- One-grain volume and free volume distributions 

Let us first concentrate on the non-vibrated smooth grains. The distribution of the 
Voronoi cell areas defined above is presented on figure |2Ja). This distribution displays 
two peaks centred on < v s >= 1.00 and < v\ >= 1.49, that corresponds respectively 




Figure 2. (a) Distribution of the Voronoi cells area for the non- vibrated smooth 
grains. Vertical dashed lines : minimal Voronoi cell area. Vertical dash-dotted 
lines : conditional average Voronoi cell area. Inset: distributions of the free volume 
conditioned by the grain size; (dark): small grains; (grey): large grains. (b) 
Distribution of the free volume for the non-vibrated smooth grains including both 
sizes of grains. Inset: density of states of one grain surrounded by a free volume v* . 
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Figure 3. Free volume distributions for the three different experimental runs: non- 
vibrated smooth (o) and rough (x) grains, and vibrated smooth grains (•) (the vertical 
lines indicate the mean values). 



to the average area occupied by the small and large grains, computed independently. 
Also indicated on the figure are the minimal values that a Voronoi cell can possibly 
take - the closest regular hexagon- for each type of grain, v™ m = \/3/2 ~ 0.866 and 
v min _ y^/2(di/d s ) 2 ~ 1.35. Both peaks present a well-defined exponential tail, which 
is easily isolated when considering the distributions of the free volume (v* gl = v — v™{ n ), 
for each type of grain as shown on the inset of figure Efa). It shall be noticed that the 
characteristic free volume for each type of grain, v{ : = 0.055 and v{* = 0.060 are closer 
to each other than one would have expected given the particle size ratio. In other words, 
the characteristic free volume accessible to one grain not only depends on the size of 
that grain, but also on its neighborhood, most often composed of grains of both sizes. 

Figure Efb) displays the distribution of the free volume without conditioning it 
to the type of grains. One still observes the well-defined exponential tail, with a 
characteristic free volume = 0.058. In the spirit of Edwards' description, one 
may assimilate this volume, characteristic of the sample, to a compactivity. This is 
of course a rather strong interpretation of the data, since it would be necessary at least 
to check that the dependence of the distributions on this compactivity is indeed fully 
embedded in the exponential tail. Still, as an indication, we have plotted in the inset 
of figure I2{b) the quantity ln(p(i/)) = ln(P(v^))) + /v? , the density of states of one 
grain surrounded by a free volume v' . It saturates for large free volumes, after a sharp 
increase for small values. 

We now characterize the dependence of the above free volume distribution on both 
the roughness of the grains and the history of the packing. Figure El displays the free 
volume distributions for the non-vibrated packings of smooth and rough grains, as well 
as for vibrated and non-vibrated packings of smooth grains. The values of the mean and 
characteristic free volumes are reported in Table ^ Obviously, the mean free volume 
decreases when vibrating the pile. It is also lower in a pile of smooth grains than in 
a pile of rough grains prepared the same way: roughness allows more efficient arch 
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building, hence larger interstitial space. The characteristic free volume associated with 
the exponential tail is smaller for the vibrated, more compact system, as suggested 
in ^j. Inside non- vibrated packings, this characteristic free volume is not much 
sensitive to the roughness of the grains. 

4-2. Volume and free volume distributions for clusters of neighboring grains 

As explained in the introduction, our main goal in the present study is to test 
experimentally the coherence of a statistical construction for granular media. To do 
so, one has to study not only one-grain quantities, but also the statistical properties 
of a subsystem made of a sufficiently large number of grains. In addition, it is clear 
that a richer information may be gained from the distributions of global quantities 
characterizing the system, rather than from simple averages of these quantities. Indeed, 
distributions of global quantities have proved to be a useful tool in the study of non- 
equilibrium systems, both from the experimental [231221 and from the numerical jTTTH ITTT] 
point of view. Moreover, it should be noticed that in granular media, one may not 
expect the N-body probability distribution for the volumes of a set of N grains to 
factorize as a product of one-grain distributions, due to possible correlations between 
grains. Accordingly, some interesting and non trivial properties may emerge from the 
statistical properties of systems consisting of iV ^> 1 grains, so that there is a strong 
motivation to determine such statistical properties experimentally. 

To this aim, we first define the system on which measurements are to be performed. 
As the packing is disordered, there is no obvious way to determine the system by a simple 
geometric boundary, contrary to what may be done if grains were sitting on a regular 
lattice. To overcome this difficulty, we define clusters of N grains by starting from an 
arbitrary grain, and looking for its N — 1 nearest neighbors, that is, the N — 1 grains 
that are the closest to the original grain. The cluster is then composed of the original 
grain and of the N — 1 nearest neighbors. 

Now that the clusters of grains are defined, one would like to characterize their 
statistical properties. At this stage, it is important to note that the free volume and the 
full volume are no longer equivalent, since clusters are composed of both small and large 
grain. Figure Ufa) and (b) respectively display the distribution of the volume per grain 
vn = A^ _1 S^ 1 t> i , and of the free volume per grain v^ N = A^~ 1 S^ 1 f/, inside clusters of 
N neighboring grains. 

The distribution of the free volume per grain v N converges much faster towards a 
well defined distribution than that of the volume. Indeed, it is clear from figure Ufa), in 
the case of the volume distributions, as iV increases from 1 to 10, that combinatory effects 
induced by the bidispersity give rise to series of peaks which complexify the distribution. 
For N = 2, one observes three peaks, which are related to the three two-particle cluster 
configurations. Similarly, for N = 3 -not shown here-, the distribution exhibits four 
peaks, etc. Such effects are absent in the case of the free volume distributions (figure^)). 
Also, describing the volume distribution imposes to know the minimal value of the 
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volume per grain, a non trivial quantity for clusters composed of grains of both sizes, 
whereas by construction, the distribution of the free volume readily addresses this issue. 
Accordingly, we restrict the following analysis to the free volume distributions. 

Given that the free volume per grain v^ N may take arbitrary positive values, and 
that the distribution P{v^ N ) seems to have an exponential tail, at least for small N, it is 
natural to try to describe our data with Gamma laws. This description is also supported 
by numerical simulations of Poisson-Voronoi cells [22. • We thus use the following ansatz 



which involves two unknown parameters t]n and xn] Vn m &y be considered as a shape 
parameter whereas xn is a scale parameter. Note that at this stage, one shall not 
anticipate any relation between xn and the compactivity as defined by equation (J2J). 
A simple way of determining rj^ and xn is to adjust their values so that the first two 
moments of the Gamma law coincide with the experimental data. To this aim, we 
make use of the relations (v^) = t] N xn an d Var(Vjy) = i] N x%, from which one deduces 
Xn = Var(v^)/(v^) and r/jy = ( v n)/Xn- The average free volume (v N ) is plotted on 
figure a), and is seen to rapidly converge to a constant value i>4> the mean free volume 
of the packing. A more surprising behavior is obtained from the variance of v N , when 
plotted as a function of 1/N, as shown on figure Efb). If the individual free volumes 
associated with each grains were independent random variables, the variance of the total 
free volume of the cluster would scale as N, so that free volume per grain in the cluster 
would be proportional to 1/N. On the contrary, for the range of cluster sizes over which 
measures have been performed, the variance Vax(v N ) scales as 

Var(v£) « vN~ a , < a < 1 (8) 

The above large N behavior of the two first moments leads us to the following scaling 
Xn = XeN~ a and tjn = VeN a , with Xe = v/v/^ and r\ e — v 4 /v. The values of the 
different parameters are reported in Table H 
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Figure 5. Dependence on N of (v^) (a) and Var(w^) (b) for the three different 
experimental runs (same symbols as in Fig.^J). Curves in (b) are shifted for clarity. 



Altogether, for large N one can rewrite the distributions P{v^ N ) in the following 
form, making the scaling with N explicit 



/ N 
P{Vn) = Z c ( Xe , Ve ,N) 



cxp 



Xe 



(9) 



(10) 



where the normalization factor Z c (x e , Ve, N) is given by 

Z c = exp N a r] e (in — - 1 
L V Ve 

Hence a non-extensive scaling appears over the explored range of cluster 
sizes for lnP(v^). Figure El displays the scaling function, &{vfj) = 
N~ a \ia[P(v^) exp(N a v^/xe)]- The collapse and the linearity of the data confirm the 
validity of the above analysis, and yields for the distribution $^(1^) of the free volume 
Vf = Nv{j for a subsystem of iV grains : 



Z c (Xe,Ve,N) 



exp 
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0.120 


0.154 
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0.059±0.001 


0.046±0.001 


0.061±0.001 


Xe 


0.048±0.0015 


0.032±0.0015 


0.040±0.0015 


Ve 


3.0±0.1 


3.7±0.1 


3.9±0.1 


a 


0.85±0.02 


0.80±0.02 


0.77±0.02 



Table 1. Summary of the different parameters extracted from the distributions of free 
volume in the three systems (same symbols as in Fig. 0). 

In the extensive case (a = 1), one would recover the general form of equation (JSJ, 
hence a formal identification of \ e with the compactivity X. Here this is not the case 
and one does not have a compactivity independent of the system size, when applying 
the definition of equation @. However, one can still define generalized entropy and 
"free energy" densities s(w^) and f{v^ N ) by 

s(v f N ) = r] e \nv f N and f(v f N ) = v f N - Xes{v f N ) (12) 

which allows P{v N ) to be rewritten as 

P(„/) = ^ e -^/(^)/xe (13) 

a form similar to the equilibrium one, except for the exponent a ^ 1. With such 
definitions one recovers, from a saddle-point calculation on equation Q, the relation 
X^ 1 — 9s/dv*\/ v f\, which generalizes the one usually used in equilibrium extensive 
statistical mechanics and suggest that Xe can still be interpreted as the compactivity of 
the packing. 

It should be noticed that the above canonical analysis implicitely assumes that the 
total volume of the packing is constant, so that, when considering a cluster of N grains, 
the rest of the packing plays the role of a reservoir. In the present finite system we 
can estimate the total volume fluctuations to be between 1.7% and 0.8% depending on 
whether we assume the sub-extensive scaling to hold up to the system size or not. This 
is consistent with a variation of half a layer of grain at the top of the pile. 

5. Discussion 

5.1. Physical origin of the non-extensivity 

As explained above, our central result is the observation of a sub-extensive scaling, with 
the cluster size N, of the variance of the free volume per grain in clusters of increasing 
sizes. This unusual scaling qualitatively indicates the presence of correlations between 
the individual free volumes of the grains, as otherwise Var(t>^) would scale as 1/N. 

Let us try to discuss more quantitatively how the departure from the 1/iV-behavior 
for Var(iriy) relates to the presence of correlations in the system. To this aim, we make 
the following simplifying assumptions: (i) the cluster has the shape of a disk, with a 
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radius R; (ii) the density of grains is homogeneous inside the cluster; (iii) introducing 
5v{ = v{ — v-f, where is the average free volume per grain, the correlation (SvfSvj) 
depends only on the relative distance r = |r$ — r 3 -|, and not on further information 
like the position of the grains % and j with respect to the center of the cluster. This 
correlation is denoted as C(r) in what follows. Expanding the variance, one finds 

Var(^) = Var ( ^£>/ ) = ]p Z>*W> ( 14 ) 
\ i / ij 

Under the above assumptions, one can then write 

Var UfEM = (^) 2 / dr ' I * 2 C(| ri -r 2 |) (15) 

\ i J J\ti\<R J\r 2 \<R 

Let us now assume that the correlation function C(r) decays like r~ 7 for large r. Then 
the large R behavior of the integral given in Eq. (fT3j) can be computed, and turns out to 
be Var(v^) ~ iV~ 7 / 2 . To obtain this result, one may use polar coordinates to integrate 
over the vectors r x and r 2 . The angular part of the integral actually appears only in the 
prefactor of the power law iV~ 7 / 2 . This suggests that the precise shape of the cluster 
only plays a minor role, and does not affect the scaling behavior itself. In the present 
case, given the values of a reported above, one would infer that correlations decay like 
r~ 7 with 7 w 1.6 over the experimentally accessible range of distances. 

From a physical point of view, such correlations seem to be related to anomalous 
fluctuations of the numbers of particles inside the clusters. More precisely, one may think 
either about anomalous fluctuations of the density of particles [IB], or to fluctuations 
of the fraction of small and large particles inside a cluster. Figure Ufa) displays the 
averaged number of grains (N) and the averaged number of small grains (N s ) inside a 
disk of radius R (note that this is a procedure different from the one previously used, 
which instead involved clusters of N neighboring grains). As expected both scale like 
R 2 confirming the visual observation that the granular packing is compact. Conversely, 
Var(iV) and Var(A^ s ), shown on figure 0(b), do not scale like R 2 , and hence not like N 
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or N s respectively. Note that such anomalous fluctuations are much stronger for the 
fraction of small particles than for the density itself. The fact that the average value 
scale precisely like R 2 is interesting, as it shows that trivial finite-size corrections due 
to, e.g., the discrete character of granular media, are very small in the range of sizes 
explored. This confirms that the anomalous scaling of the variance is relevant, and is 
not simply due to generic finite-size corrections affecting any observable. 

Such anomalous fluctuations and scalings could be originated in the two- 
dimensional confinement of our system, and one may wonder whether the same exponent 
would be expected in three dimensions. Such experiments have not been conducted in 
three dimensions. However one element of answer goes as follows. We ran simulations 
of Poissonian Voronoi tesselations in two dimensions (that is, with point-like particles 
thrown at random on a plane), performed the same size scaling analysis and did not 
observe any anomalous scaling. This underlines the role of the finite extension of the 
grains in the observed scalings. 

5.2. Interpretation of the parameters rj e and Xe 

Interestingly, even if extensivity does not hold, the above experimental study still leads 
to a consistent statistical description. In particular, the relation x~ l = ds/dv*\/ v f\, 
which generalizes the one usually used in equilibrium extensive statistical mechanics, 
may be recovered from a saddle-point calculation on equation This results indeed 
gives support to the interpretation of \ e as the compactivity of the packing. It is also 
consistent with the definition proposed for thermodynamic parameters like temperature 
in equilibrium systems with long range interactions |SS] • In this case, these parameters 
are essentially rescaled (with respect to more usual definitions) by a size- dependent 
effective number of degrees of freedom, so that they remain finite in the thermodynamic 
limit. We also note that the qualitative behavior of Xe is the one physically expected. 
For instance, when lowering the volume fraction by vibrating the packing, Xe is seen to 
decrease, in a way similar to what happens in an equilibrium system, where temperature 
decreases when lowering energy. 

The second parameter describing the statistics of the packing, r] e , rather 
characterizes the entropy s(vl). In more physical terms, f] e N a may be thought of as the 
effective number of independent (collective) objects in a cluster of N grains. This may 
be seen from the most simple case of the Poisson- Voronoi cells in one dimension [32J. 
In this case, points are thrown at random on a line, and the length of the associated 
Voronoi cells is given by half the sum of the length of the two segments surrounding 
the point considered. As the lengths of the segments between two adjacent points are 
independent random variables, a Voronoi cell may be thought of as having two internal 
degrees of freedom, thus suggesting r\ = 2. This simple reasoning is confirmed by exact 
computations [32]. So coming back to our experimental situation, this suggests that 
the presence of correlations leads to complex and effectively independent objects, with 
possibly a fractal shape, the number of which is growing non-extensively with the size 
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of the cluster. 

5. 3. Consequences for the thermodynamic formalism 

The above results have important consequences from the point of view of measurements 
of statistical properties in granular systems. On the one hand, they show that it is 
important to check the extensivity of the system when extracting the compactivity from 
usual thermodynamic relations. The non-extensive scaling behavior is expected to hold 
up to a scale comparable to the correlation length of the packing. It has been recently 
proposed that amorphous media may be in a state close to a critical state associated 
with the "jamming transition" |2*2*1 123| 124*]. Hence, this suggests that the correlation 
length may be quite large in disordered granular packings. 

Other consequences also need to be stressed. The parameter rj e appearing in the 
expression of the entropy varies together with Xe when the way to prepare the system 
is changed. This is a signature of the non-equilibrium character of the system. In an 
equilibrium situation, changing the preparation conditions would amount to change, say, 
the energy and the temperature of the system, but it would not affect the functional 
expression of the entropy. Only the value of the entropy would change. Here the value of 
rj e is determined by the microscopic dynamics and the change in r\ e suggests that the way 
microscopic configurations are sampled by the preparation process changes drastically 
when varying the density: not only the typical volumes sampled are modified, but also 
the way configurations of a given volume are explored. 

It is worth underlining that such a situation could formally occur in the case of 
Edwards' approach. Indeed, since Edwards' equiprobability assumption is restricted 
to the states blocked under the considered dynamics, the density of such states may a 
priori depend on that dynamics. From that point of view, the generalized formalism 
proposed in Sect. 2, takes into account any possible departure from the equiprobability 
of microstates, including the specific case of a dynamical selection of the blocked states. 
The drawback is then that the probability weights f(C) are unknown, but the relevant 
macroscopic information may be encoded into the parameter r] e , which is experimentally 
accessible. 

Finally, we note that the presence of correlations "renormalizes" in some way the 
"compactivity" Xe, so that it is no longer equal to the characteristic volume v** of the 
exponential decay of the one-grain distribution. Although one expects that extensivity 
is recovered for system sizes larger than a crossover length scale related to the correlation 
length, there is no reason to believe that this "large scale compactivity" would be equal 
to v'*. Hence, this shows that measuring the slope of the exponential tail of the one- 
grain distribution is not a relevant determination of compactivity. 
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6. Conclusions 

In this paper, we have studied clusters of increasing sizes in a bidimensional granular 
packing. Our central result is that a sub-extensive scaling of the fluctuations of 
free volume appears. This strongly suggests the presence of correlations between the 
individual free volumes of distant grains, over the range of sizes explored experimentally. 
In spite of this nonextensive scaling which partly rules out the statistical formalism 
proposed for granular materials, it is still possible to define in a consistent way an 
intensive compactivity. This parameter is found to differ significantly from the one 
that may be inferred from the exponential tail of the one-grain distribution. In some 
sense, this leads to the idea that correlations act as a "renormalization" of the statistical 
parameters. ^From a more general point of view, it also shows that one has to be careful 
when applying statistical concepts to a system for which the microscopic dynamics is 
poorly understood -namely a kind of "black box" . In order to be sure that the procedure 
is well enough controlled, one should make some checks of the fundamental properties 
of the system, among which the extensivity plays a major role. 
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